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Abstract 

We study the problem of all sky search in reference to continuous gravitational 
wave (CGW) whose wave-form are known in advance. We employ the concept of Fit- 
ting Factor and study the variation in the bank of search templates with different Earth 
azimuth at t = 0. We found that the number of search templates varies significantly. 
Hence, accordingly, the computational demand for the search may be reduced up to 
two orders by time shifting the data. Keywords: gravitational wave, data analysis, all 
sky search, LIGO. 
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1 Introduction 



The ground based laser interferometric gravitational wave (GW) detectors viz. LIGOHl, 
VIRGOQ, GEO600I3 1 TAMA 300||1 produce a single data stream that may contain con- 
tinuous, chirp, burst and stochastic GW. These detectors don't point, but rather sweep their 
broad quadrupolar beam pattern across the sky as the earth moves. Hence the data analysis 
system will have to carry out all sky searches for its sources. In this, the search of contin- 
uous gravitational wave (CGW) without a priori knowledge appears to be computationally 
quite demanding even by the standard computers expected to be available in the near fu- 
ture. It appears that due to limited computational resource it will be not feasible to do all 
sky all frequency search for the CGW in the months/year data set. Hence for the search, 
one have to look for the signal in the short duration of data set. Understanding the problem 
and the target sensitivity of the advance LIGO, it may be feasible to do all sky search for 
one day data set at low and in narrow frequency band. The search may be more significant 
if it is done in the frequency band where most of the Pulsars are detected by other means. 
Also, the choice of optimal data processing and a clever programming is also integral part 
of a solution to this problem. Amongst these the pre-correction of time series due to the 
Doppler modulation before the data is processed may be a method, which will reduce the 
computational requirements. In reference to this, Schutz||5l has introduced the concept of 
patch in the sky as the region of space throughout which the required Doppler correction 
remains the same. He shown that the number of patches required for 10^ sec. observation 
data set of one KHz signal would be about 1.3 x 10^^. This also implies that the bank 
of search templates require for the match filtering in an all sky search. However, the size 
of the patch would also depend on the data analysis technique being employed, which in 
turn depend on the parameters contain in the phase of the modulated signal. Hence, in this 
paper after incorporating the azimuth of the Earth (Port) at t = in the Fourier transform 
(FT) obtained by Srivastava and Sahay|6|, we investigate its effect in the bank of search 
templates for one sidereal data set. Hence in the next section we extend the FT of the Fre- 
quency modulated (FM) CGW obtained by themQ with taking account of Port- In section 
3 we employ the concept of Fitting Factor (FF)|I3 and check the cross correlation of the 
templates with the corresponding data set exceeds the preassigned threshold by considering 
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the source location as parameters at different (3orb of the signal manifold and compute the 
number of search templates require for an all sky search apphcable to such analysis. We 
present our conclusion in section 4. 

2 Fourier transform of the Frequency modulated contin- 
uous gravitational wave 

The FT analysis of the Frequency modulated CGW has been done by Srivastava and 
Sahay[6 1 by taking account the effects arising due to the rotational as well as orbital motion 
of the Earth. However, they have neglected an important parameter Port- To obtain FT by 
taking account of f3orb, we rewrite the phase of the received CGW signal of frequency fo at 
time t with some modification given by them|l6| and may be written as 



$(t) = 271 fot + Z COs{a^rot - (t) +N'cOs{^rot - 5) - M (1) 

where 



(2) 



M = ^ {RseSm9cosa + VWTWcos6) , 

Z = ^RseSinO, AT = ^y/WTW, 

V = -ReSinQ;(sin^sin0cose + cos6'sine) , 

Q = Re sin a sin 9 cos (j) , 

a = - /Sorb , S = tan"^ ^ - (3rot , 

a = Worb/Wrot ~ 1/365.26, Worbt = a^rot, 

n = (sin 9 cos 0, sin 6' sin , cos 9) , ^rot = Wrott 



where 9, 0, R^., Rse, Worb Wrou a and e represent respectively the celestial co-latitude, lon- 
gitude. Earth radius, average distance between Earth centre from the origin of SSB frame. 
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orbital and rotational angular velocity of the Earth, co-latitude of the detector and obliquity 
of the ecliptic. Here (3orb and (3rot are the azimuth of the Earth and detector at t = respec- 
tively. 

Now, the two polarisation functions of CGW can be described as 

= Vcos[$(t)] (3) 

h^{t)=ho,smmt)] (4) 
where /io+, ho^ are constant amplitude of the two polarizations. 

Considering the function 

h{t) = cos[<l>(t)] (5) 
the FT for one sidereal day may be given as 



h{f) = I "'^ cos[<l'(t)]e-^^^^*dt ; Tabs = one sidereal day = 86164 s (6) 







d 

which may be split into 

kf)] =/.-+4+; (7) 

- d 

2Wrot Jo 

— / g-i[5v++^cos«-(T)+A^cos(C-5)-A^]^^ ^g-j 
/o ± / 

= ; i = ^rot = Wrott (10) 

Jrot 

As I^^ contributes very little to h{f) . Hence, hereafter, we drop from Eq. ^ and 

L J d 

write u in place of . Using the identity 



l=oo 

giiKcos^ = J^{±k) +2j2t^Ji{±K)cosW (11) 

1=1 
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we obtain 

Hf) , 



iM 



2tt 



k=oc 



k=l 



2Wrot 

JoiAf) + 2 ^ JmW"" cosm(e - 6) 



X 



Jo{Z) + Jk{Z)i^ cos k{a^ - a) 



m=l 



where J stands for the Bessel function of the first kind. After integration we get 

k=oo m=oo 



Hf) 



V 



d 2w< 



rot 



k=—oo m=—oo 



(12) 



(13) 



A 
B 
C 
V 



(k+m)iT 



u^ — {ak+m)^ 

sin2z/7r cos(2aA;7r — ka — m5) — |cos 2z/7r sin(2afc7r — ka — m5) + sm{ka + m5)} 
cos 2z/7r cos(2aA;7r — ka — m5) + gj^ 2utx sin(2aA;7r — ka — m5) — cos{ka + m5) 



Now its straight forward to obtain the FT of the two polarisation states of the wave and can 
be written as 



h+{f) = KAh{f) 



fc=oo m=oo 



2Wr: 



rot 



J2 E ^'Mc-^^ 



(14) 



fc=— oo m=— oo 



-^ho. \kf) 



k=oo m=oo 



2w. 



rot 



(15) 



fc=— oo m=— oo 



For the analysis we reduce the computing time ~ 50% by using symmetrical property of 
the Bessel functions by rewriting the Eq. (1131 as 
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Hf) 



V 



Wrot 



^^^^^^^4^[{sm-M -sin(A^ -2z/7r)} + z{cos - cos(A^ - 2z/7r)}] + 



"=°° 7 (\f) 

Jo{z) T [{yu - XV) - i{xu + yv)] + 



m=l 
k=oo m=oo 

k=l m=—oo 



(16) 



X = sm{M-mn/2) 

y = cos{M-mn/2) 

U = sin2z/7rcosm(27r — 5) — ^ {cos2z/7rsinm(2z/7r — 5) — sinm(5} 

V = cos 2z/7r cosm(27r — 5) + — sin 2z/7r sinm(27r — 5) — cosm5 



(17) 



3 Bank of search templates for an all sky search 

The study of templates has been made by many research workers ll8l l9l fTOl fTT i . However, 
the question of possible minimum efficient interpolated representation of the correlators for 
an all sky search is a problem of interest. In this, the study of the variation in the bank of 
search templates for the short duration of data set in reference to the parameter contain in 
the phase of modulated signal is very important. Hence we check the effect of Porb in the 
modulated signal by plotting h{f) (Figures ([T]) and ©) for LIGO detector at Hanford 

L J d 

(the position and orientation of the detector can be found in Ref. 12) of unit amplitude 
signal for 



fo = 50Hz, Po = 0,tt/2 
9 = n/18, = 7r/4 



(18) 



Here we take the ranges lfT3lfT4l of k and m as 1 to 27300 and -10 to 10 respectively 
because the value of Bessel functions decreases rapidly as its order exceed the argument. 
From the plot, we observe the obvious but major shift in the spectrum with Port- Under- 
standing the effect we check the variation in the bank of search templates at different Porb- 
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The bank of search templates are discrete set of signals from among the continuum of 
possible signals. Consequently all the signals will not get detected with equal efficiency. 
However, it is possible to choose judiciously the set of templates so that all the signals 
of a given amplitude are detected with a given minimum detection loss. FF is one of the 
standard measure for deciding what class of wave form is good enough and quantitatively 
describes the closeness of the true signals to the template manifold in terms of the reduction 
of SNR arising due to the cross correlation of a signal outside the manifold with the best 
matching templates lying inside the manifold, given as 



where h{f) and hxif] 0^, (px) represent respectively the FTs of the actual signal wave form 
and the templates. The inner product of two waveform hi and /12 is defined as 



where * denotes complex conjugation, ~denotes the FT of the quantity underneath. S'„(/) 
is the spectral noise density of the detector and has been taken stationary and Gaussian as 
the bandwidth of the signal is extremely narrow. 

To compute the number of search templates we consider the LIGO detector at Han- 
ford, receive a CGW signal of frequency fo = 50 Hz from a source located at (9, 0) = 
(0.1°, 30°). We chosen the data set such that (3orb = at t = 0. In this case we take the 
ranges of k and m as 1 to 310 and —10 to 10 respectively and bandwidth equal to 2.0 x 10"^ 
Hz for the integration. Now, we select the spacing A^^ = 4.5 x 10^^, thereafter we max- 
imize over by introducing spacing A0 in the so obtained bank of search templates and 
determine the resulting FF. In similar manner we obtain the FF at Porb = and n/2. 
The results obtained are shown in the Fig. ©• We observe that the nature of the curve is 
similar. Hence, it may be interesting to obtain a best fit of the graphs and may be given as 

Nremplates = 10^^ [Cq + CiX - C2X^ + CsX^ - C^X^ + C^X^ - CqX^ + Cyx'^] ; (21) 



FF 



{h{f)\hT{f;9TAT)) 



(19) 



VihAf; 0T, 0t)|/^t(/; ^t, <l>T)){Hf)\hif)) 





(20) 
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0.80 <x< 0.995 
where Co, Ci, C2, C3, C4, C5, Ce, Cy are the constants given in the Table ©. 

In view of the above investigation, the grid spacing (A6', A0) in the (6*, 0)— parameter 
of templates may be expressed as 

A^^ = T{FF, U e, 0, T„,„ f3,rb) (22) 
Similarly, 

A0 = g{FF, /„ e, 0, T,fc„ (23) 

In Ref. to 15, for one sidereal data set, the dependence of FF on the template variables 
9t and (pT is given as 

pp = g-0.00788(e - St)" (24) 



= g-O.O1778(0 - 4>T? (25) 

From Eqs. dH), <E1I) and (El), we may write 

J^(FF,50,0.r,30°,lrf,/3o,b) = [-(0.00788)-Mn(FF)]^^^ (26) 



^(FF, 50, O.r, 30°, Id, Porb) = [-(0.01778)-^ In(FF)] (27) 

Hence, for the selected FF one can determine A9 and Acf). However, there is no unique 
choice for it. Here we are interested in the assignment of A^^ and A(f) such that the spacing 
is maximum resulting into the least number of templates. As we have mentioned earlier, 
there is stringent requirement on reducing computer time. Accordingly, there is serious 
need of adopting some procedure/formalism to achieve this. For example, one may adopt 
the method of hierarchical search. 
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Table 1: Coefficients of tlie best fit graphs obtained for the bank of search templates. 



(3orb Co Cl C2 C3 C4 C5 Cq C7 



0° -4.36537 34.4523 116.44 218.462 245.734 165.719 62.0404 9.94642 



45° -403.012 3187.44 10795.5 20296.6 22877.7 15459.9 5799.45 931.641 



90° -263.622 2086.48 7071.80 13305.6 15008.9 10150.3 3810.65 612.645 



4 Conclusions 

We have incorporated the parameter Port in the FM signal and investigated its effect in the 
bank of search templates for an all sky search. The analysis for complete response of the 
detector has not been done as the requisite is analogous to what presented by them|6| and 
also not required for the templates analysis. 

We observe that the change of Port affects the spectra severely. Consequently in the 
bank search templates. On investigation for one sidereal day data set of a signal of 50 Hz, 
we found that the number of search templates varies significantly with Port- For the case 
investigated we found that for FF = 0.97 approximately 24.8267 x 10^°, 22.7840 x 10^^ 
and 32.3097 x 10^^ search templates may be require, if Port = 0, 7r/4 and n/2 respectively. 
Hence, the computational demand may be reduce up to two orders by time shifting the 
data. One may optimize the require number of search templates by doing analysis with the 
method given by Owen lfT6l . 
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Figure 1 : Spectrum of a Doppler modulated signal when (3orb — 0. 




Figure 2: Spectrum of a Doppler modulated signal when Po^b = 7r/2. 
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Figure 3: Variation in the number of search templates with FF at at different /3, 
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